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This  eontribution  eoneerns  a  general  purpose  fluid-multibody  system  (MBS)  simulation  framework  that  ean 
be  used  to  analyze  the  fluid-solid,  two-way  eoupled  dynamies  at  low  to  medium  Reynolds  numbers  (0  <  Re  < 
1500).  The  simulation  framework  ean  be  leveraged  to  investigate  MBS  applieations  that  inelude  (/)  rigid  and 
flexible  bodies  of  arbitrary  geometry;  (//)  bilateral  eonstraints  (joints);  (///)  unilateral  eonstraints  assoeiated  with 
impaet  and  eontaet  phenomena;  and  (/v)  frietion/eohesion. 

The  fluid  dynamies  problem  is  formulated  using  the  fluid  momentum  and  eontinuity,  i.e.,  Navier-Stokes 
equations.  These  equations  are  spatially  diseretized  via  a  weakly  eompressible  smoothed  partiele  hydrodynamies 
(SPH)  Lagrangian  method  [1],  whieh  relies  on  moving  markers  to  store  state  information  assoeiated  with  fluid 
phase.  The  spaee  dependent  variables,  sueh  as  veloeity  and  pressure,  are  smoothed  out  loeally  via  a  sealar  fiinetion. 
That  is,  to  obtain  a  variable,  a  gradient,  or  a  hydrodynamies  foree  at  an  arbitrary  loeation  of  the  domain,  one  needs 
to  aeeount  for  partial  eontributions  eoming  from  nearby  markers.  External  forees  sueh  as  fluid-solid  interaetion 
(FSI)  foree  are  added  to  the  hydrodynamies  foree.  The  fluid  equations  of  motion,  whieh  upon  spatial  diseretization 
beeome  a  set  of  ordinary  differential  equations,  are  solved  explieitly  using  a  seeond  order  Runge-Kutta  integration 
method. 

Of  several  approaehes  that  have  been  eonsidered  in  the  literature  to  model  the  fluid-solid  eoupling,  we  show 
that  using  a  point-eloud  diseretization  of  a  solid  results  in  an  aeeurate  ealeulation  of  the  fluid-solid  eoupling  forees 
[2].  In  this  approaeh,  the  MBS  dynamies  is  solved  by  providing  the  solver  with  distributed  forees  eaptured  by  the 
point  eloud  representation.  Alternatively,  a  eonstraint -based  approaeh  was  examined  where  the  non-penetration 
eonditions  were  treated  as  unilateral  eonstraints  on  fluid  markers.  In  the  latter  approaeh,  the  rigid  body  forees  were 
not  marker  resolution  independent. 

The  solution  strategy  outlined  has  been  implemented  in  Chrono  as  a  dedieated  add-on  ealled  Chrono::FSI  [3]. 
Figure  1  shows  a  vehiele  model  used  in  a  fording  simulation  along  with  its  point  eloud  representation.  The 
refinement  of  the  point  eloud  ean  be  based  on  the  aeeuraey  required  in  the  fluid  system.  The  vehiele  is  modeled 
in  Chrono: :  Vehiele,  whieh  is  a  Chrono  toolkit  that  provides  a  template  based  approaeh  to  rapid  vehiele  prototyping. 
In  Chrono:: Vehiele,  templates  define  the  basie  modeling  elements  (bodies,  joints,  foree  elements),  impose  the 
subsystem  topology,  preseribe  the  design  parameters,  and  implement  the  eommon  funetionality  for  a  given  type 
of  subsystem  (e.g.  suspension)  partieular  to  a  speeifie  template  (e.g.  double  wishbone). 

A  heterogeneous  parallel  eomputing  approaeh  is  adopted  to  speed  up  the  simulation.  The  state  of  the  fluid 
system  is  stored  and  updated  on  a  graphies  proeessing  unit  (GPU),  while  the  MBS  dynamies  solution  relies  on  a 
shared-memory,  multi-eore  implementation.  The  effieieney  gains  obtained  by  adopting  a  heterogeneous 
eomputing  approaeh  ean  be  hindered  by  the  massive  data  transfer  between  the  GPU  and  host,  whieh  takes  plaee 
via  a  relatively  high  lateney,  low  bandwidth,  PCIe  bus.  We  eireumvent  this  issue  by  minimizing  the  data  transfer 
to  and  from  the  GPU.  To  this  end,  the  point  eloud  data  for  eaeh  rigid  body  resides  in  the  GPU  memory.  After 
ealeulating  the  distributed  forees  on  the  point  eloud,  those  forees  are  redueed  in  parallel  using  the  thrust  library 
[4].  The  resulting  forees  and  torques  (six  variables  per  rigid  body),  are  transferred  to  the  MBS.  Similarly,  after 
updating  the  MBS,  the  GPU  updates  the  loeations  and  veloeities  of  the  points  of  the  eloud,  relying  on  the  new 
position,  orientation,  veloeity,  and  angular  veloeity  of  the  rigid  objeets. 

Chrono:  :FSI  has  been  used  for  vehiele  mobility  in  fording  operations  as  shown  in  Figure  2.  The  eomputational 
time  per  simulation  time  step  is  almost  entirely  dietated  by  the  fluid  eomponent  and  is  similar  to  that  of  a 
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simulation  of  a  pure  fluid  system  of  the  same  size.  This  indieates  that  the  minimal  data  transfer  between  the  GPU 
and  the  host  system  does  not  slow  down  the  simulation.  Finally,  it  is  worth  mentioning  that  Chrono::FSI  allows 
for  the  simulation  of  a  broad  elass  of  fluid-MBS  interaetion  problems  that,  beyond  the  fording  analysis  diseussed, 
eoneern  the  general  purpose  dynamies  of  rigid  and/or  flexible  bodies  interaeting  with  fluid. 


Fig.  1:  The  vehicle  model  used  in  a  fording  simulation.  Any  component  that  interacts  with  the  fluid,  such  as  the  vehicle,  tire  geometries, 
and  the  solid  boundaries  shown  in  the  left  image,  are  represented  as  a  point  cloud  in  the  fluid  system.  The  resolution  of  the  point  cloud 
discretization  can  be  changed  based  on  the  accuracy  or  the  specific  dynamics  required. 


Fig.  2:  Simulation  of  the  vehicle  mobility  in  fording  using  Chrono::FSI 
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